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ABSTRACT 

Herbig-Haro (HH) jets are commonly thought of as homogeneous beams of 
plasma traveling at hypersonic velocities. Structure within jet beams is often 
attributed to periodic or "pulsed" variations of conditions at the jet source. Sim- 
ulations based on this scenario result in knots extending across the jet diameter. 
Observations and recent high energy density laboratory experiments shed new 
light on structures below this scale and indicate they may be important for under- 
standing the fundamentals of jet dynamics. In this paper we offer an alternative 
to "pulsed" models of protostellar jets. Using direct numerical simulations we ex- 
plore the possibility that jets are chains of sub-radial clumps propagating through 
a moving inter-clump medium. Our models explore an idealization of this sce- 
nario by injecting small (r < rj et ), dense (p > p je t) spheres embedded in an 
otherwise smooth inter-clump jet flow. The spheres are initialized with velocities 
differing from the jet velocity by ~ 15%. We find the consequences of shifting 
from homogeneous to heterogeneous flows are significant as clumps interact with 
each other and with the inter-clump medium in a variety of ways. Structures 
which mimic what is expected from pulsed-jet models can form, as can previ- 
ously unseen "sub-radial" behaviors including backward facing bow shocks and 
off-axis working surfaces. While these small-scale structures have not been seen 
before in simulation studies, they are found in high resolution jet observations. 
We discuss implications of our simulations for the interpretation of protostellar 
jets with regard to characterization of knots by a "lifetime" or "velocity his- 
tory" approach as well as linking observed structures with central engines which 
produce the jets. 

Subject headings: ISM: jets and outflows - ISM: Herbig-Haro objects - hydrody- 
namics 
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1. Introduction 

Herbig-Haro objects have been the subject of significant analytical, observational, and 
numerical attention since their discovery. Observations using optical (e.g. Bally et al. 2002) 
and IR (e.g. Velusamy et al. 2007) techniques reveal that these jets typically show striking 
large scale collimation extending out to parsec distances combined with features appearing 
on a range of smaller scales. 

Structure along the jet beam ("knots" or "clumps") have, in particular, received consid- 
erable attention. The origin of knots remains a subject of debate. Early studies focused on 
clumpyness of the HH bow shocks; Norman & Silk (1979) postulated the existence of single 
"interstellar bullets," while Schwartz (1978) attributed the structures to stationary clumps 
being overrun by a wind. Stationary crossing shocks due to an overpressured jet beam ex- 
panding and then re-collimating were an early possibility that was considered for knots along 
the beam (Biihrke et al. 1988; Raga et al. 1990a). More recently, Rubini et al. (2007) have 
suggested oblique shock focusing as a natural mechanism for hydrodynamic knot formation, 
though the presence of magnetic fields (Hartigan et al. 2007), precession (Masciadri et al. 
2002), and interactions with the environment (Raga et al. 2002a; de Gouveia dal Pino 1999; 
Yirak et al. 2008) all offer other means by which dense clumps might be created. While 
considerable work has gone into these scenarios, currently the most favored model for the 
knots are internal working surfaces where shocks are driven down the beam by pulsation at 
the jet source. This "pulsation" model was first proposed by Rees (1978) and was exten- 
sively explored by Raga and collaborators (Raga et al. 1990b; Biro & Raga 1994; Raga & 
Biro 1993). In pulsed jet simulations the density and velocity cross-sectional profiles Pj{r) 
& Vj(r) in the jet-launching region are kept fixed, while the magnitude of the velocity varies 
sinusoidally (Raga et al. 1990b; Volker et al. 1999). The pulsation scenario has become so 
dominant that even when attempting to address questions unrelated to clump formation, 
periodic inflow variations are frequently employed (e.g. Suttner et al. 1997). 

A variety of observational signatures can be recovered via pulsed jet models through 
careful choice of specific jet physical parameters and sinusoidal variability. In Raga et al. 
(2002b), for example, a two-mode launching model was proposed using velocity histories ex- 
tracted from observations of HH 34 and HH 111. Using these pulsation modes axisymmetric 
hydrodynamic simulations provided a convincing match to the location of the leading bow 
shock and the location of bright knots in the beam. These and similar results provide strong 
support for pulsed jet models. 

A detailed examination of jets observed at the highest spatial resolution however shows 
features which do not fit into the pulsed jet paradigm. In particular a number of "archetypal" 
jets show features at scales below the jet radius (r < rj) which are distinctly displaced from 
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the jet axis. In the case of HH 47 the jet clearly shows a non-axisymmetric morphology 
in the form an apparent helical bending of the beam (Hartigan et al. 2005). The beam 
itself is defined by a sequence of quasi-periodic knots with displacements to either side of 
the nominal jet axis. Explanations for this bending have included impacts with objects 
(Raga et al. 2002a; de Gouveia dal Pino 1999), magnetic fields (Hartigan et al. 2007) and 
precession of the jet source (Masciadri et al. 2002). We consider the presence of sub-radial, 
non-axisymmetric features to be a challenge to the pulsed jet paradigm. Here, we consider 
an alternative to the pulsed jet model and investigate the consequences of intrinsic density 
heterogeneity in jet formation and evolution. 

We are motivated to explore this model both by observations and new high energy den- 
sity laboratory astrophysics (HEDLA) plasma experiments. Using pulsed power wire array 
technologies Lebedev et al. (2005), Ciardi et al. (2007), and Ciardi et al. (2008, in prepara- 
tion) have presented experiments that track the evolution of fully magnetized, hypersonic, 
radiative jets. The stability of hydro and MHD jets has long been a topic of debate, and 
these experiments shed some light on the real dynamics of 3-D systems (Xu et al. 2000). 
The experiments show that kink mode instabilities strongly affect the jet. As the kink mode 
grows into the non-linear regime its disrupts but does not destroy the jet. The saturation of 
the instability transforms the jet into a sequence of collimated chains of knots which prop- 
agate with a range of velocities. Similar fragmented chains have been seen in other pulsed 
power experiments (Golingo et al. 2005). 

Thus we propose a model in which the velocity and density profile of the jet are variable 
in time and space on scales less than the jet age and jet radius, respectively. Specifying 
conditions on these scales allows the model to achieve complex structures not seen before in 
simulations. Our simulations utilize AMR techniques and so are able resolve clumps in the 
jet beam at acceptable levels. 

In the following sections, we describe the model in detail, offer dynamical and observa- 
tional signatures, and briefly discuss the results. 

2. Computational Method & Physical Model 

Numerical simulations were undertaken with the AstroBEAR computational code 1 . 
AstroBEAR is a parallel adaptive mesh refinement (AMR) code which allows a variety of 
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choices for numerical solvers, integration schemes, and cooling modules for hydrodynamic 
or magnetohydrodynamic astrophysical fluids (Cunningham et al. 2006, 2007). Here, the 
code solved the 3D hyperbolic system of equations for inviscid, compressible flow using a 
spatial second-order and temporal first-order accurate MUSCL scheme using a Roe-averaged 
linearized Riemann solver. Simple radiative cooling is included separately using an iterative 
source term with the cooling curve of Dalgarno & Mccray (1972). 

In what follows, the geometry of the system is taken to have the jet axis aligned with the 
z-axis, with x- and y-axes following the right handed convention. A base grid of 36 x 36 x 96 
cells covered the simulation domain of extent 1, 200 x 1, 200 x 3, 200 AU. AstroBEAR employs 
a patch-based adaptive grid to refine areas of interest; here, two levels of refinement were 
used, yielding a maximum effective resolution of 288 x 288 x 768 cells. This corresponds 
to a minimum cell length of Ax = 6.23 x 10 13 cm (~4 AU). All boundaries had outflow 
conditions, with user-specified conditions in the jet launching region. 

The ambient number density and temperature were initialized to constant values of 
p a =10 cur 3 and T a =2,000 K, respectively, giving a sound speed of c a = 5.26 km s -1 . The 
inhomogeneities in the beam were introduced as spherical density and velocity perturbations 
in an otherwise smooth beam: we shall refer to the former as "clumps" and the latter 
as the "jet." The jet had number density, radius, velocity, and temperature of pj=10 2 
cm -3 , rj=100 AU, Vj—150 km s _1 , and Tj=2,000 K. This resulted in a Mach 30, over-dense 
(Xja = Pj/Pa = 10), and overpressured (pj/p a = 10) jet. The jet was launched at x,y=600 
AU on the z = plane. The clumps all had the same initial number density of p c =10 3 cm" 3 , 
yielding density ratios of Xcj = Pc/Pj=10 and Xca = Pc/Pa=10 2 . The clumps were seeded 
with radii and velocities that were random within ranges of r c =26-60 AU, and t> c =132-168 
km s -1 . Relative to the jet, these correspond to ranges r c /r.,=0.26-0.60 and t> c /t>,,=0.88-i.l2. 
The x- and ^/-locations of the clumps in the jet were random with the constraint that the 
entire clump be located within the jet beam. All clumps were seeded at the same z location 
in the grid. The clumps had temperatures such that they were in pressure balance with the 
jet beam (assuming a relative velocity of zero). They were seeded nominally 8 years apart; 
this resulted in the production of 12 clumps before the end of the simulation (t sim = 100 yr). 

At the maximum AMR level, the jet radius was resolved by 24 cells, and the clump radii 
by 6-14 cells. This represents probably the lower limit on desirable resolution. Jet, clump, 
and ambient materials were separately tracked with passive advected tracers. 
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3. Results & Analysis 



A time-sequence of the simulation is given over four panels in Fig. 1. Density plots and 
a Schlieren image are shown in Fig. 2. Figure 1 shows a 3-D representation of the simulation 
in the form of a set of iso-density contours. In the panels, the jet beam enters from the left 
hand side of the grid and propagates to the right. Shortly after the start of the simulation, 
knots appear with random sizes, locations, and speeds. The figure has been adjusted to track 
the evolution of the clumps via an iso-density contour of a passive clump tracer (in green). 
Thus the clumps are readily recognizable as initially spherical inclusions within the beam 
close to the inflow boundary cells. As the simulation progresses, the clumps evolve via their 
interaction with the inter-clump material in the beam and, in some cases, with other clumps. 
We note that bow shocks form for clumps with high enough differential speeds relative to 



Fig. 1. — Isocontours of logarithmic density at four times in the simulation, t =30, 53, 77, & 
100 yr. The clumps are depicted in light green, with the jet material in blue. The x-z plane 
along the jet axis clips the jet material contours. 

the beam i.e. |At> c | = \v c — Vj\ > Cj where Cj is the sound speed in the jet beam. As a 
bow shock propagates into the beam material, a transmitted shock will propagate into and 
compress the clumps. If the differential velocity 8v c = Av c / Vj is positive then the bow shock 
faces forward into the direction of jet propagation. If Sv c is negative then the bow shock will 
face backwards towards the jet source. These bow shocks are not directly visible in Fig. 1, 
however the sign of 8v c is apparent in the compression occurring on the leading (trailing) 
edges when Sv c is positive (negative). The bow shocks are apparent in Fig. 2. In the bottom 




-6 - 



panel, the clump at z — 2100 AU has 5v c > and hence a forward facing shock forms at its 
leading edge while the larger clump at z — 900 AU shows the opposite behavior. We note 
also the presence "spur shocks" in the bottom panel of Fig. 2. These are concave shocks that 
originate at the edges of the beam and arc away from the beam (Heathcote et al. 1996). We 
find the presence of clumps off-axis naturally leads to the development of such structures. 

Once clumps are launched into the jet there are three possible consequences. The first 
consequence is the clump propagating downstream unimpeded and colliding with the jet 
head. In this case the dense clump may break through the bow shock defining the front 
edge of the jet, leading to significant non-axisymmetric structures there. This behavior is 
apparent in both Fig.'s 1 and 2 where a dense clump has already traversed the jet length 
and propagated through the jet shock/bow shock structure at the terminus of the beam. 
The presence of a significant "knob" protruding at the lower edge of the jet head defines the 
extent of the clump which now forms the leading edge of the jet. 

A second possibility however is that the clump will not make it to the leading edge 
of the jet. The behavior of shocked clumps has been extensively studied both analytically 
and numerically (e.g. Klein et al. 1994; Jones et al. 1996; Poludnenko et al. 2002; Fragile 
et al. 2004). These studies demonstrate that a clump compressed by a strong transmitted 
shock wave (in this case the transmitted wave originates from its relative motion within 
the beam) will eventually be destroyed in a "cloud crushing time" given approximately by 

1 /2 

t cc = 2r c x C j /|Ai> c | where Xcj — Pel Py If this happens before the clump reaches the jet 
head then the clump material will be dispersed within the beam. Thus the cloud crushing 
timescale should be compared with the timescale required to cross the length of the jet beam, 
tj C = Lj/(v c — Vf, s ). Note that Lj = v^t where v\, 8 is the speed of the jet head given by the 
familiar formula Vf, s = Vj/(1 + Xja^) an( f t is the time at which the clump is launched. 
Comparing these expressions we find the critical launching time t*, 



A clump of radius r c , velocity dispersion Av c and density ratio Xcj needs to be launched at 
a time t = ti aunc h < t* in order for it to reach the jet head before being destroyed. For each 
instantiated clump, comparing ti aunc h to t* reveals that t\ aunc h < t* for three of the twelve 
clumps. This implies that late in time (assuming no other interactions) nine of the twelve 
clumps would disperse before reaching the jet head. However, t cc > t S i m for all the clumps 
in the simulation (t sim is the duration of the simulation). We therefore expect all clumps 
to exert a strong influence on the jet beam throughout the simulation, and we expect some 
of the clumps to affect strongly the morphology of the jet head. These expectations are 
confirmed in the simulation. 





(1) 
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The third possibility for the long term evolution of a clump, assuming 5v c is not the 
same for all clumps, is interaction with another clump. The interaction can take the form 
of direct or glancing collision depending on the impact parameter b. Even when b > 2r c 
there can still be interactions between a clump-driven bow shock and a neighboring clump 
(Poludnenko et al. 2002) Such interactions will be determined by the width of the clump bow 
shock which will, in general be determined by the Mach number of the clump through the 
inter-clump beam media (M c = \Av c \/cj). The collisions of clumps is also a process which 
has been well studied, and one expects the formation of transmitted shock waves within 
the clumps which heat and compress clump material as it streams into the shock (Klein 
et al. 1994; Miniati et al. 1999). Figure 1 and Fig. 2 show a number of such interactions 
occurring. By the last panel of Fig. 1, clump collisions have resulted in a merged structure 
near the head of the jet, and their effect on each other and on the jet beam itself is complex. 
It is noteworthy however that the collision, compression and subsequent merger of clumps 
can come to resemble the internal working surfaces in homogeneous pulsed jets. Note the 
resemblance in the bottom panel of Fig. 2 of the structure at z = 2300 AU to the thin shock 
bounded working surfaces seen in jet pulsation simulations. An important difference in the 
case of clumps however is that these structures remain sub-radial. The merged structure 
identified at z = 2300 AU has a lateral width of h w 1.8rj and is offset from the jet axis by 
tq 0.25rj. Thus, even collisions between slow moving clumps overtaken by faster moving 
ones differ from what is expected for homogeneous jet beams with varying inflow velocity. 
In this case, the resulting shock structures are slightly smaller than the jet diameter and 
displaced from the jet axis. 

The difference between homogeneity and heterogeneity is particularly striking for glanc- 
ing clump collisions: r c < b < 2r c . In these cases the clump-clump interaction will be off 
center and one can expect from momentum conservation that non- axial motions will result. 
The top two panels of Fig. 2 illustrate this point, showing the off-center collision of three 
clumps. Before the collision the velocity vectors of all three clumps is purely axial v = v z e z . 
After the collision the clumps have acquired transverse v r velocities. The ability to generate 
non-axial motions within the beam via clump interactions is an important point as proper 
motion studies of highly resolved HH jets show knot to knot variability in both direction 
and speed (Hartigan et al. 2005; Bally et al. 2002). Moreover, this interaction forms the 
structure at z — 2300 AU seen in the bottom panel of Fig. 2 which therefore has varying 
velocity components across its surface, in contrast to similar structures in pulsed-jet mod- 
els. We note again that the clump injection time, position within beam cross-section and 
velocity were all chosen randomly within constraints. We now address the issue of the ob- 
servational properties of jets with this kind of heterogeneity, with respect to assumptions of 
pulsation. In many studies lifetimes of features along an HH beam jet beam are derived by 




1000 
800 




Fig. 2. — The top two panels give grayscale images of logarithmic density in the x-z plane 
located on the jet axis at two different times, t =53 & 77 yr, corresponding to the lower 
left and upper right panels of Fig. 1. Lighter gray corresponds to denser material. Velocity 
vectors originating at three knots have been overlaid, and they are seen to change as the 
knots interact. The bottom panel shows a synthetic Schlieren image at t — 100 yr, which 
illuminates such features as a clump with a forward bow shock at z =2100 AU, a clump with 
reverse bow shock at z =900 AU, and clump-induced "spur shocks" at several places along 
the jet beam. The disk-like feature at z — 2300 AU is discussed further in the text. 
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relating the current position of the feature and its proper motion (Bally et al. 2002; Hartigan 
et al. 2005). Going a step further several some studies (e.g. Raga & Noriega-Crespo 1998; 
Raga et al. 2002b)) infer full velocity histories of jets from observations in terms of multi- 
ple pulsation modes at the jet launch region. While these studies are able to convincingly 
reproduce some observational characteristics of the jets, our models of heterogeneous jets 
shed new light on the issue of recovering pulsation histories from observations. We note 
that inferring ejection histories from current positions and velocities requires two intrinsic 
assumptions: first, that the launching is smoothly varying and periodic in nature, and sec- 
ond, that pulsation-formed knots are coherent throughout their lifetimes. In particular, it 
is assumed that given a position zq for a hypersonic, semi-ballistic jet, the time-variability 
in the past can be inferred from the current observed velocity structure. Considering the 
current epoch to be t = one uses observed axial velocity v (z ) to compute dynamical times 
to = — zq/v(zq) for fluid parcels when they where ejected from the source. The velocity 
history of the jet as Vj(t) — v [t — —z/v(z)} can then be derived. In Raga et al. (2002b) for 
HH 34 and HH 111 it was assumed that Vj(t) could be described by a periodic function of 
the form 



HH 111, n = 1, while for HH 34, accounting for the collimated knots near the source 
requires n — 3. 

Assuming Eq. 2 for Vj(t) can bias the description of jets as we now demonstrate. For 
each data frame in our simulation, we may perform a similar analysis; see Fig. 3. Along the 
axis of the jet we may denote the locations of "knots" as those positions where the density 
peaks at some value above the jet density. The axial density profile is given in the left panel 
of Fig. 3, with knots identified by circles. We then take the velocities at these positions to 
be the knot velocities. This allows us to calculate dynamical times and, using a single mode 
(n = 1) for simplicity, reconstruct a separate velocity history for each data frame in the 
simulation. As implied by the second assumption above, in order for this analysis to be valid 
the resulting amplitude and period should be robust to the time in the simulation, especially 
late in time when there are the greatest number of knots. We found that the results were 
not robust, as depicted in the right hand panels of Fig. 3. The three panels show, from top 
to bottom, the result of least-squares fits to the data at t — 100, 80, and 40 yr with 1-a 
error bounds overlaid. The top panel shows the fit resulting from the corresponding data 
given in the left panel. The average velocity of the knots remains roughly constant for each 
of the three times depicted and is close to the jet velocity (vj = 150 km s -1 ), as expected. 
However, the amplitude varies by ~ 30% over time, and the period varies by over 100%. 





where Vj , 



r k and (f>k are the amplitude, period, and phase of each pulsation mode. For 
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Fig. 3. — Left panel: The axial density profile normalized by the ambient density, p/p a , 
is plotted with positions designated as "knots" given by "o". Right panels: From top to 
bottom, velocities for knots ("°") versus their launching (dynamical) time at t =100, 80, & 
40 yr. Results from least-squares fits are shown with l-a error overlaid. The period r and 
mean velocity with single-mode sinusoidal amplitude, Vj = Vj$ ± vf\ are printed on each 
panel. 
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The goodness of the fits, though less at t = 80 yr than the other two times depicted, appears 
adequate. Although not explicitly shown here, variation in the fit results shows correlation 
with knot-knot interactions. 

4. Discussion and Conclusions 

We have performed 3-D simulations using a new "collimated clumps" scenario for pro- 
tostellar jets. Our model offers a fundamentally different paradigm for understanding jet 
origins and dynamics in the sense that heterogeneity is seen as being intrinsic and links 
the jet morphology on "meso-scales" to the processes (such as instabilities) occurring on 
"micro-scales" near the central engine. 

While the pulsed-jet model has been successful at interpreting some aspects of jets, it 
may be misleading if used too generally. In particular, the assumption of sinusoidal pulsa- 
tions can limit the interpretation of HH object observations. In our simulations, which had 
no sinusoidal variation in time, we nonetheless were able to recover (erroneous) sinusoidal 
behavior using an analysis similar to that which has been carried out in the past on obser- 
vations. While one part of this behavior (the mean velocity Vj$) fit the initial conditions 
well, the assumption of periodic pulsation allows the false conclusion that the structures in 
the beam arose due to periodic ejection behavior. We conclude therefore that care should 
be taken when attributing observations of apparent sinusoidal velocity in protostellar jets to 
corresponding sinusoidal behavior of a central engine. 

In contrast to pulsed-jet models, our model offers two attractive features: first, a natural 
mechanism (knot-knot interactions) helps explain small-scale features along the jet axis. 
The idea that knots or bow shocks in HH objects are evolving clumpy structures has been 
discussed before in the context of observations (e.g. Reipurth et al. 2002), with ongoing 
HST observations giving continued support to this idea (e.g. Hartigan et al. 2005). The 
interaction of distinct sub-radial knots with each other and with an overall bow shock offers 
a simple explanation for such evolution. The second feature the model offers is the presence 
of unique observational characteristics in the form of forward- and reverse-facing bow shocks 
and spur shocks at the edges of the jet beam. Our scenario provides a simple mechanism 
for the formation of sub-radial non-axisymetric features via multiple dense, clumps which 
have a nonzero velocity dispersion. The ease with which non-axisymmetric features like spur 
shocks (Heathcote et al. 1996) develop in our models is attractive. 

What would be the origin of the entrained clumps? The issues of the stability of jets 
particularly at "micro-scale" regions near the launch region remains unresolved (Konigl 2004; 
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Ouyed et al. 2003; Xu et al. 2000; Micono et al. 1998). In addition, issues of magnetic 
field strength and velocity perturbations near the jet launch region also present problems 
(Hartigan et al. 2007). Although we do not explicitly consider details of disk/jet launching 
here, plausible scenarios present themselves. It has long been considered that the interaction 
between the forming star and accretion disk — via magnetic fields — is responsible for jet 
launching (Konigl 1982). Further, anisotropy of the magnetized disk during the YSO phase is 
a distinct possibility, as suggested by Combet & Ferreira (2008). One could therefore theorize 
the occurrence of a non-cylindrically-symmetric accretion burst as a possible explanation for 
non-axial density enhancements in the jet beam. 

Recent high energy density laboratory astrophysics (HEDLA) investigations provide 
an unique window into the behavior of fully 3-D radiative hypersonic MHD jets. These 
experiments demonstrate that magnetized jet beams in the lab may break up into a sequence 
of quasi-periodic knots due to the kink instability (Ciardi et al. 2007). These knots may 
be displaced slightly to the side of the nominal jet axis and may propagate with varying 
velocities. This results in morphologies qualitatively reminiscent of HH-jet beams. It should 
be noted that the present simulation does not employ magnetic fields; however, it remains 
an open question whether magnetic fields remain dynamically important on the length scales 
of consideration here (Hartigan et al. 2007; Ostriker et al. 2001). It seems plausible that a 
similar process could occur in the astrophysical context, beginning with a smooth beam near 
the central engine which then becomes disrupted owing to the kink or sausage instabilities on 
intermediate scales. This would result in a series of knots which continue to evolve as they 
propagate away from the central engine. Such a scenario would also explain the observed 
velocity differences between knots, attributable to the particulars of each knot's formation. 
The present simulation is an idealization of this model. 

The degree to which the results of the simulation presented here depend on the values 
chosen for parameters such as Xcj and Av c requires further investigation. In particular, 
tiaunch < t* was not satisfied for all of the instantiated clumps, implying that some should 
disperse before reaching the jet head. The present simulation did not progress far enough 
to witness this behavior (that is, t cc > t sim ). Also, using a Av c ^ and varying from 
clump to clump should result in an injection of considerable vorticity into the jet beam. 
Future work should therefore explore the long term evolution of the clumped jet. The aspect 
ratio (length/width) of the jet's final state in the present simulation is less than that for HH 
objects, implying again that simulations which progress for longer times would be of interest. 
While increased resolution would of course be of benefit, extracting details of clump-clump 
interaction within the jet-beam environment is probably outside the scope of this study 
and would be appropriate for a separate investigation (Dennis et al. 2008, in preparation). 
Finally, the inclusion of magnetic fields to study scenarios of knot formation and the resulting 
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observational consequences would be of additional interest. 

We would like to thank Pat Hartigan, Sergey Lebedev, and Andrea Ciardi for their time. 
Tim Dennis and Brandon Shroyer also provided invaluable support and help. 

Support for this work was in part provided by by NASA through awards issued by 
JPL/Caltech through Spitzer program 20269 and 051080-001, the National Science Founda- 
tion through grants AST-0507519 as well as the Space Telescope Science Institute through 
grants HST-AR-10972, HST-AR-11250, HST-AR-11252 to. We also thank the University of 
Rochester Laboratory for Laser Energetics and funds received through the DOE Cooperative 
Agreement No. DE-FC03-02NA00057. 

REFERENCES 

Bally, J., Heathcote, S., Reipurth, B., Morse, J., Hartigan, P., & Schwartz, R. 2002, AJ, 123, 
2627 

Biro, S., & Raga, A. C. 1994, ApJ, 434, 221 

Biihrke, T., Mundt, R., & Ray, T. P. 1988, A&A, 200, 99 

Ciardi, A., Ampleford, D. J., Lebedev, S. V., & Stehle, C. 2008, ApJ, 678, 968 

Ciardi, A., Lebedev, S. V., Frank, A., Blackman, E. C, Chittenden, J. P., Jennings, C. J., 
Ampleford, D. J., Bland, S. N., Bott, S. C, Rapley, J., Hall, G. N., Suzuki- Vidal, 
F. A., Marocchino, A., Lery, T., & Stehle, C. 2007, Phys. Plas., 14, 6501 

Combet, C, & Ferreira, J. 2008, A&A, 479, 481 

Cunningham, A. J., Frank, A., & Blackman, E. G. 2006, ApJ, 646, 1059 

Cunningham, A. J., Frank, A., Varniere, P., Mitran, S., & Jones, T. W. 2007, ArXiv:astro- 
ph/0710.0424, 710 

Dalgarno, A., & Mccray, R. A. 1972, ARA&A, 10, 375 

de Gouveia dal Pino, E.M. 1999, ApJ, 526, 862 

Fragile, C. P., Murray, S. D., Anninos, P., & van Breugel, W. 2004, ApJ, 604, 74 
Golingo, R. P., Shumlak, U., & Nelson, B. A. 2005, Phys. Plas., 12 



- 14 - 

Hartigan, P., Frank, A., Varniere, P., & Blackmail, E. G. 2007, ApJ, 661, 910 

Hartigan, P., Heathcote, S., Morse, J. A., Reipurth, B., & Bally, J. 2005, AJ, 130, 2197 

Heathcote, S., Morse, J. A., Hartigan, P., Reipurth, B., Schwartz, R. D., Bally, J., & Stone, 
J. M. 1996, AJ, 112, 1141 

Jones, T. W., Ryu, D., & Tregillis, I. L. 1996, ApJ, 473, 365 

Klein, R. I., Mckee, C. F., & Colella, P. 1994, ApJ, 420, 213 

K6nigl, A. 1982, ApJ, 261, 115 

Kdnigl, A. 2004, ApJ, 617, 1267 

Lebedev, S. V., Ciardi, A., Ampleford, D. J., Bland, S. N., Bott, S. C, Chittenden, J. P., 
Hall, G. N., Rapley, J., Jennings, C. A., Frank, A., Blackman, E. G., & Lery, T. 2005, 
MNRAS, 361, 97 

Masciadri, E., de Gouveia dal Pino, E.M., Raga, A. C, & Noriega-Crespo, A. 2002, ApJ, 
580, 950 

Micono, M., Massaglia, S., Bodo, G., Rossi, P., & Ferrari, A. 1998, A&A, 333, 1001 
Miniati, F., Jones, T. W., & Ryu, D. 1999, ApJ, 517, 242 
Norman, C, & Silk, J. 1979, ApJ, 228, 197 

Ostriker, E. C, Stone, J. M., & Gammie, C. F. 2001, ApJ, 546, 980 

Ouyed, R, Clarke, D. A., & Pudritz, R. E. 2003, ApJ, 582, 292 

Poludnenko, A. Y., Frank, A., & Blackman, E. G. 2002, ApJ, 576, 832 

Raga, A., & Noriega-Crespo, A. 1998, AJ, 116, 2943 

Raga, A. C, Binette, L., & Canto, J. 1990a, ApJ, 360, 612 

Raga, A. C, Binette, L., Canto, J., & Calvet, N. 1990b, ApJ, 364, 601 

Raga, A. C, & Biro, S. 1993, MNRAS, 264, 758 

Raga, A. C, de Gouveia dal Pino, E.M., Noriega-Crespo, A., Mininni, P. D., & Velazquez, 
P. F. 2002a, A&A, 392, 267 

Raga, A. C, Velazquez, P. F., Canto, J., & Masciadri, E. 2002b, A&A, 395, 647 



- 15 - 

Rees, M. J. 1978, MNRAS, 184, 61P 

Reipurth, B., Heathcote, S., Morse, J., Hartigan, P., & Bally, J. 2002, AJ, 123, 362 
Rubini, F., Lorusso, S., Del Zanna, L., & Bacciotti, F. 2007, A&A, 472, 855 
Schwartz, R. D. 1978, ApJ, 223, 884 

Suttner, G., Smith, M. D., Yorke, H. W., & Zinnecker, H. 1997, A&A, 318, 595 

Velusamy, T., Langer, W. D., & Marsh, K. A. 2007, ApJ, 668, L159 

V61ker, R., Smith, M. D., Suttner, G., & Yorke, H. W. 1999, A&A, 343, 953 

Xu, J., Hardee, P. E., & Stone, J. M. 2000, ApJ, 543, 161 

Yirak, K., Frank, A., Cunningham, A., & Mitran, S. 2008, ApJ, 672, 996 



This preprint was prepared with the AAS IAT^X macros v5.2. 



